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60/267,804, filed February' 8, 2001, which is hereby fully incorporated by reference 
herein as though set forth in full. 

1. Field of the Invention. 

This invention relates generally to the field of barycentric coordinates for 
1 5 representing interior points of polygons in terms of their boundary vertices, and, more 
specifically, methods of computing such coordinates, and applications of such 
coordinates. 

2. Related Art. 

Barycentric coefficients are useful for representing interior points of a polygon 
20 in terms of the polygon vertices. Barycentric coordinates are weights which, when 
applied to the vertices of a polygon, and the results added, yield the interior point. 

Methods for computing triangular barycentric coordinates have been known 
for centuries, but efforts to generalize these methods to n-sided polygons have not 
been entirely successful, in part because the formulations resulting from these efforts 
25 are subject to one or more weaknesses which limit their effectiveness and 
applicability. 

For example, in Eugene L. Wachspress, "A Rational Finite Element Basis", 
Mathematics in Science and Engineering, Volume 1 14, Academic Press, 1975, a 
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method for computing barycentric coordinates for n-sided polygons is disclosed 
which can be explained with reference to Figure 1. 

Referring to Figure 1, n-sided polygon 100 is illustrated, and it is desired to 
compute the barycentric coordinates for the interior point p, identified in the figure 
5 with numeral 102, in terms of the vertices of polygon 1 00. As is known, for an n- 
sided polygon, there will be n barycentric coordinates for an interior point, one 
corresponding to each of the n vertices. These barycentric coordinates, which can be 
referred to as a , V/, j e [l...n\ are such that the following relationship holds: 

f:j 10 where aj refers to the barycentric coordinate applied to the vertex <jj. (As used herein, 
a bolded item refers to a vector, while an unbolded item refers to a scalar). 

According to the technique proposed by Wachpress, a weight Wj for vertex qj 
if; is derived from the product of the area of the triangle 106 (shaded in Figure 1) formed 
5 " by the three adjacent vertices qj+i, qj, and qj_i and the areas of the n-2 interior 
Ci 15 triangles, identified with numerals 108, 1 10, and 1 12, formed by the point p and the 
polygon's adjacent vertices (but excluding the two interior triangles that contain the 
vertex qj). Mathematically, this formulation may be expressed as follows: 

w, =4q j _ 1 ,q j ,q j+1 ).n ts!{/ , y+l} 4q k -i,q I c 5 p) (2) 
where the function A refers to the area of the triangle whose vertices are the 
20 arguments of the function. 

This formulation is then repeated for each of the vertices in the polygon. The 
resulting weights may then be normalized in accordance with the following 
expression to form the barycentric coordinates for the point p: 

(3) 



This formulation, however, is complex since the computation of the 
barycentric coordinate corresponding to a particular vertex depends on all the vertices 
of the polygon, and entails repetitive computations of triangle areas. Consequently, 
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the formulation is computationally intensive and therefore unacceptable for 
applications where a faster approach is needed. 

Summary 

In a first embodiment, the invention provides a method for deriving 
barycentric coordinates for a point p within an n-sided polygon. For a particular 
coordinate Wj, corresponding to the vertex <jj, this method embodies a formula which 
depends only on the edge pqj, and the two adjacent angles 8j and y } . This formula 
may be expressed as follows: 

~ Mil 2 

This computation may be repeated for each of the n polygon vertices to result in 
barycentric coordinates corresponding to each of the n vertices. If necessary, the 
coordinates may then be normalized using expression (3). 

In one implementation, for accuracy reasons, the cotangent function used in 
this formula is not computed through a trigonometric function call. Instead, in this 
implementation, the cotangent is derived from a division between the dot product of 
two vectors derived from the triangle involved, and the cross product of these two 
vectors. In particular, for a triangle having vertices abc, the cotangent of the angle at 
vertex b may be derived using the following formulation: 

bc * ba (5) 

||bcxba| K) 

where ba is the vector a - b, and be is the vector c - b. 

In a second embodiment, the invention provides a method for deriving weights 
Wy for expressing a vertex i in a mesh representation of an object surface in terms of 
its one-ring neighbors j e N{i). For a particular vertex qy, this method embodies a 
formula which depends only on the edge qiqj, and the two adjacent angles 8j and yj. 
This formula may be expressed as follows: 

= cot(rJ+cotfo) (6) 
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This computation may be repeated for each of the vertices in the mesh representation 
to result in weights allowing each of the vertices to be expressed in terms of its one- 
ring neighbors. Normalized weights may then be derived from non-normalized 
weights using an expression analogous to expression (3). 

In a third embodiment, the invention provides a method of parameterizing a 
mesh representation of an object surface using the formulation expressed in (6). This 
method begins with the step of computing the weights wy in 3D space (in contrast to 
parameter space) for each of the vertices in the mesh representation. For a vertex i, 
i e the weights w y allow the vertex i to be expressed in terms of its one-ring 

neighbors j e N(i) . 

Once the weights have been derived, the method proceeds to fixing the 
positions of boundary vertices in parameter space, such as by assigning each of these 
vertices a position on a fixed boundary C. For a boundary vertex i, the position 
assigned to this vertex in parameter space may be referred to as C u . 

Then, the method proceeds to solving the following system of linear equations 
in order to derive the parameterization of the mesh representation: 

Wy (uj - Uj ) = 0 if i is an interior vertex 

jeNO) \ (7) 

Uj = C u if i is a boundary vertex 

where u f is the vertex i in parameter space (and uj is the vertex j in parameter space), 
and C u is the boundary position in parameter space assigned to the boundary vertex i. 

The result of solving this system of linear equations is the coordinates Uj of 
each of the vertices i in parameter space, / e [l ...n] . 

Other systems, methods, features and advantages of the invention will be or 
will become apparent to one with skill in the art upon examination of the following 
figures and detailed description. It is intended that all such additional systems, 
methods, features and advantages be included within this description, be within the 
scope of the invention, and be protected by the accompanying claims. 



Vi,ie[l. ...n\ 
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Brief Description of the Drawings 

The components in the figures are not necessarily to scale, emphasis instead 
being placed upon illustrating the principles of the invention. In the figures, like 
5 reference numerals designate corresponding parts throughout the different views. 

Figure 1 illustrates a representation of an n-sided polygon partitioned into 
triangles useful for computing the barycentric coordinates of a point p using the prior 
art formula of Wachpress. 

Figure 2 illustrates a representation of an n-sided polygon showing the edge 
10 pqj and the adjacent angles 8j and yj useful for computing the barycentric coordinates 
of a point p using a method according to the invention. 

Figure 3 illustrates parameterization of a mesh representation of an object 
surface. 

Figure 4 illustrates application of a method according to the invention to 
1 5 computing weights for expressing a vertex of a mesh representation in terms of its 
one-ring neighbors. 

Figure 5 illustrates the vectors and triangle vertices which are useful in 
computing the cotangent of one of the vertex angles according to a method which 
avoids a trigonometric function call. 
20 Figure 6 is a simplified flowchart illustrating one embodiment of a method of 

parameterizing a mesh representation of an object surface using a formula according 
to the invention. 

Detailed Description 

25 As utilized herein, terms such as "about" and "substantially" are intended to 

allow some leeway in mathematical exactness to account for tolerances that are 
acceptable in the trade, e.g., any deviation upward or downward from the value 
modified by "about" or "substantially" by any value in the range(s) from 1% to 20% 
of such value. 
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Referring to Figure 2, in a first embodiment of the invention, a method for 
deriving non-normalized barycentric coordinates for a point p within an n-sided 
polygon is provided. For a particular coordinate w j5 corresponding to the vertex qj, 
identified with numeral 202, this method embodies a formula which depends only on 
5 the edge pqj, identified with numeral 204, and the two adjacent angles 8j and yj, 

identified respectively with numerals 206 and 208. This formula may be expressed as 
follows: 



cot^)+cot(<5j) 




(8) 



This computation may be repeated for each of the n polygon vertices to result in non- 
1 o normalized barycentric coordinates corresponding to each of the n vertices. 

Normalized coordinates may then be derived from the non-normalized coordinates 
using the following expression: 




(9) 



Compared to the Wachpress formulation (2), this formulation has the 
15 advantage that it is local and hence computationally simple: only the edge pqj and its 
two adjacent angles 5j and Vj are needed. 

The coordinates derived from this formula possess several desirable 
properties. First, they satisfy the following affine combination property: 

p= 2 a j*l> with pf' =1 (10) 
20 This property may be expressed in the following form using the non-normalized 
coordinates w*: 

5>,(qj-p)=0 (ID 

Mi--] 

Second, the coordinates {a, } =1 n (and the normalized coordinates) are 
smoothly dependent on p and the vertices of the polygon. This ensures smoothness in 
25 the variation of the position of p responsive to any movement in the vertex qj. 
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Third, if the polygon is convex, the coordinates {a, } =1 satisfy the following 

convex combination property: 

^>0V/e[l...«] (12) 
This guarantees no under- and over-shooting in the coordinates: all the coordinates 
5 will be between zero and one. 

In one implementation, for accuracy reasons, the cotangent function used in 
the formula (7) is not computed through a trigonometric function call. Instead, in this 
implementation, the cotangent is derived from a division between the dot product of 
two vectors derived from the triangle involved, and the cross product of these two 
10 vectors. In particular, referring to Figure 3, for a triangle having vertices abc, the 
cotangent of the angle 0 at vertex b, which angle is identified with numeral 302, may 
be derived using the following formulation: 

be • ba 

be ba (13) 
where ba is the vector a - b, and be is the vector c - b. 
15 In one example, the formula (8) may be implemented using the following 

pseudo-code: 

// Compute the barycentric weights for a point p in an n-gon Q 
// Assumes p is strictly within Q and the vertices qj are ordered. 
computeBarycentric(vector2d p, polygon Q, int n, real w[ ]) 
20 weightSum = 0 

foreach vertex qj of Q: 

prev = (j + n - 1) mod n 
next = 0 + 1) mod n 

Wj = (cotangent(p, qj, q prer ) + (cotangent(p, qj, q^)/||p - q,| 2 
25 weightSum += Wj 

// Normalize the weights 
foreach weight w/. 

Wj 1= weightSum 

30 // Compute the cotangent of the non-degenerate triangle abc at vertex b 
cotangent(vector2d a, vector2d b, vector2d c) 
vector2d ba = a - b 
vector2d be = c - b 
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f be • ba ^ 

return ^ — T[ 

ybcxba||J 

The formula (8), any derivation thereof, and any instructions or program code 
implementing the same, may be tangibly embodied, through storage or otherwise, on 
5 or in a processor readable medium or memory including but not limited to RAM, 
ROM, PROM, EPROM, EEPROM, hard disk, floppy disk, CD-ROM, DVD, flash 
memory, etc. For purposes of this disclosure, the term "memory" refers to any 
medium which is capable of storing data. 

Referring to Figure 4, in a second embodiment of the invention, a method is 
10 provided for deriving weights wy for expressing a vertex qi in a mesh representation 
400 of an object surface in terms of its one-ring neighbors qj, V/' e N(i) . For a 
particular vertex qj, identified with numeral 402, and a neighbor vertex qj, identified 
with numeral 402, this method embodies a formula which depends only on the edge 
qiqj, identified with numeral 404, and the two adjacent angles Sj and yj, identified 
15 respectively with numerals 406 and 408. This formula may be expressed as follows: 
w = ^) + cot(^) (H) 

h-*A 2 

This computation may be repeated for each of the one-ring neighbors of q, resulting in 
weights allowing the vertex qs to be expressed in terms of its one-ring neighbors. 
Normalized weights may then be derived from non-normalized weights using an 
20 expression analogous to expression (9). 

This series of computations may be repeated for each of the other vertices in 
the mesh representation, thus allowing each vertex to be expressed in terms of its one- 
ring neighbors. 

Thus, in contrast to the previous embodiment which concerned expressing an 
25 interior point of an n-sided polygon in terms of its vertices, in this embodiment, a 
vertix of an n-sided polygon may be expressed in terms of its one-ring neighbors, 
which may be vertices of the same or a different polygon. 
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The formula (14), any derivation thereof, and any instructions or program 
code implementing the same, may be tangibly embodied, through storage or 
otherwise, on or in a processor readable medium or memory including but not limited 
to RAM, ROM, PROM, EPROM, EEPROM, hard disk, floppy disk, CD-ROM, DVD, 

5 flash memory, etc. 

In a third embodiment, the invention provides a method of parameterizing a 
mesh representation of an object surface using the formula (14). Referring to Figure 
5, parameterization is a transformation of a mesh representation 500 from 3D (x, y, z) 
space to parameter (u, v) space, and is performed to ease data manipulation and to 

10 achieve a better sampling of the original object surface than is achieved through the 
3D mesh representation. The mesh representation, after being mapped into the 
parameter space, is identified with numeral 502. 

Figure 6 illustrates a flowchart of the method. Referring to Figure 6, the 
method begins with step 602, which comprises computing the weights wy in 3D space 

1 5 (in contrast to parameter space) for each of the vertices in the mesh representation. 
For a vertex q b i e [l ...«], the weights wy are derived through application of the 
formula (14). 

Once the weights have been derived, the method proceeds to step 604, which 
comprises fixing the positions of boundary vertices in parameter space, such as by 
20 assigning each of these vertices a position on a fixed boundary C. For a boundary 
vertex i, the position assigned to this vertex in parameter space may be referred to as 
C u . Thus, in the example illustrated in Figure 5, the contour C is fixed in parameter 
space, and each of the boundary vertices of the mesh representation 500 is assigned a 
position on the contour C. For example, the vertex 504 may be assigned the position 
25 5 08 on contour C, and the vertex 506 may be assigned the position 5 1 0 on contour C. 

Then, the method proceeds to step 606, which comprises solving the following 
system of linear equations in order to derive the parameterization of the mesh 
representation: 
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Vf,/e[l....wJ 

where m is the vertex i in parameter space (and uj is the vertex j in parameter space), 
and C u is the boundary position in parameter space assigned to the boundary vertex i. 

5 The result of solving this system of linear equations is the coordinates u } of 

each of the vertices in parameter space, i e [l...w]. In the example of Figure 5, this 
step yields the coordinates of the vertex 512 in parameter space (which coordinates 
are identified with numeral 514). 

For additional information on parameterization, and the application of formula 

10 (14) to parameterization, the reader is referred to Mathieu Desbrun et al., "Intrinsic 
Parameterization of Surface Meshes," pp. 1-8, submitted for publication in 
EUROGRAPHICS 2002, and hereby fully incorporated herein by reference as 
through set forth in full. See also Mathieu Desbrun et al., "Generalizing Barycentric 
Coordinates to Irregular N-gons," pp. 1-5, submitted for publication in 

1 5 EUROGRAPHICS 2002, and hereby fully incorporated herein by reference as 
through set forth in full. 

The foregoing method may be tangibly embodied, through storage of 
instructions, program code, or the like implementing the method or otherwise, on or in 
a processor readable medium or memory including but not limited to RAM, ROM, 

20 PROM, EPROM, EEPROM, hard disk, floppy disk, CD-ROM, DVD, flash memory, 
etc. 

While various embodiments of the invention have been described, it will be 
apparent to those of ordinary skill in the art that many more embodiments and 
implementations are possible that are within the scope of this invention. 



^ w y (u j - U j ) = 0 if i is an interior vertex 
Uj = C u . if i is a boundary vertex 



(15) 



